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Recently it was found that dark energy in the form of phantom generalized Chaplygin gas may 
lead to a new form of the cosmic doomsday, the big freeze singularity. Like the big rip singularity, 
the big freeze singularity would also take place at a finite future cosmic time, but unlike the big rip 
singularity it happens for a finite scale factor. Our goal is to test if a universe filled with phantom 
generalized Chaplygin gas can conform to the data of astronomical observations. We shall see 
that if the universe is only filled with generalized phantom Chaplygin gas with equation of state 
p = — c 2 s 2 /p a with a < — 1, then such a model cannot be matched to the data of astronomical 
observations. To construct matched models one actually need to add dark matter. This procedure 
results in cosmological scenarios which do not contradict the data of astronomical observations and 
allows one to estimate how long we are now from the future big freeze doomsday. 

PACS numbers: 98.80.Cq, 04.70.-s 



I. INTRODUCTION 



Till recently there were only two types of singularities in cosmology: The big bang and the big crunch. The latter 
singularity appeared in the universe if it would go through a contracting stage which is possible both in closed universes 
(filled with baryon matter and radiation) and in open universes filled with scalar field with negative potential [l[. 

This picture changed dramatically with the coming of the dark energy revolution. Now we know about many new 
scenarios for the end of the universe among which there is the big rip singularity 0] and the sudden future singularity 
Q. Another way to obtain a cosmological doomsday was suggested in According to this scenario our vacuum 
should be rather unstable and should decay within 20 Gyr (which is possible if the gravitino is superheavy). 

A new type of future singularity, the so-called big freeze singularity, has been quite recently considered ||. It is 
a type III singularity in the notation of Ref. @ . Such a type of doomsday is possible if the universe is filled with 
phantom generalized Chaplygin gas (which was originally introduced in Q, @|, Q). Like the big rip singularity, this 
singularity would also take place at a finite future cosmic time, but unlike the big rip singularity, it happened for a 
finite scale factor. This feature results in very interesting conclusions some of which are discussed in this paper. 

The aim of this article is to find what bounds on the future lifetime of the universe filled with phantom generalized 
Chaplygin gas can be placed by current cosmological observations. As we shall see in Sec. 2 and Sec. 3, if the universe 
is only filled with generalized phantom Chaplygin gas with equation of state p = —c 2 s 2 /p a with a < —1, then such 
a model cannot match the data of astronomical observations. To construct models compatible with observations one 
need to add dark matter which allows one to obtain a class of good models where the total lifetime of the universe 
can be calculated. It is interesting to note that whereas the present age of the universe to is weakly depending on the 
parameter of the equation of state w and the Chaplygin parameter a, the total lifetime of the universe tf and the 
time to be elapsed from now to the final big freeze singularity (tf — to) are very sensitive to these parameters. For 
example, if a = —2 then for the flat universe with w = —1.01 we have tf — to = 10.3 Gyr, to = 14-4 Gyr, whereas for 
the case with w = —1.05 we have tf — to = 6.4 Gyr and t = 13.8 Gyr. 
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II. THE MODEL 



Let us consider a universe filled with dark matter (p = 0) and a generalized Chaplygin gas with equation of state 
p = —c 2 s 2 p 1+€ / 3 = — c 2 /3 2 a^p 1+£ / 3 . The the Friedmann equations would read 

a 6 J a z 
a 4ttG / ip c 2 A 

where we have denoted p = /? -6 / e (a^ — a e ) -3 / c and p — — c 2 @ 2 a e fp 1+e / 3 for, respectively, the energy density and 
pressure of the generalized Chaplygin gas. 
Observations demand that 



where k ~ 0.7. Using (|1I2I3[1 we get 



— — k a Hl, (3) 
ao 



However, p — ~ c 2 a^/3 6 / c (ay- — Aq) 1 3 / e , therefore we can find the unknown constant /3 given by 



(5) 



where Nq = -L. The current parameter of the equation of state is Wq = -f 2 - , so that 



wo = = "Z 3 Noa-aPa = -— 7 (6) 



P0 n 2 AT 3 e/3 ^0 

Finally, we have for the current vacuum energy density 



iV - 1 [kc 2 

Po = 

Using flT} and we can get c 2 ,: 



and 



_ po = N -I f kc 2 
Pc ~ ZN Q \a 2 H 2 

Using © we can derive from © the second formula for Q,q 



^^- = -^[^2+^ + 1). (9) 



3w \at H o 



We introduce now the parameter Qd,o — It is obvious that 



^d,o = 1 + 77 — 2 ~~ 



3 



Because wq = — jjrj < the condition <? d > will hold for 

fcc 2 + ff 2 fl 2 (2fc + l) 
W0< 3 (*:<»» + JT*a§) 

Let us put ko — 0.7 and consider the value of derived from (fTPl) for various values of wq and \ c „^ . The results 
attained by such a calculation are given in the following tables. 
1. k = 



w 


ft 


-1.01 


0.792 


-1.02 


0.784 


-1.03 


0.777 


-1.04 


0.769 


-1.05 


0.762 


-1.1 


0.727 



2. k = +1 



w 


Slo for -ffo a o = ^ 


fio for -ffo a o = 5c 2 


fto for H'^a z = lQc^ 


-1.01 


1.122 


0.858 


0.825 


-1.02 


1.111 


0.850 


0.817 


-1.03 


1.100 


0.841 


0.809 


-1.04 


1.090 


0.833 


0.801 


-1.05 


1.079 


0.825 


0.794 


-1.1 


1.030 


0.788 


0.758 



3. k = -1 



Wq 


fl for -ffo a o — ^ 


n Q for H'^a^ = 5c 2 


f2 for H'^a^ = 10c 2 


-1.01 


0.462 


0.726 


0.759 


-1.02 


0.458 


0.719 


0.752 


-1.03 


0.453 


0.712 


0.744 


-1.04 


0.449 


0.705 


0.737 


-1.05 


0.444 


0.698 


0.730 


-1.1 


0.424 


0.667 


0.697 



III. INTEGRATION PROCEDURE 



We have 



d 2 8ttG I /T 6 / £ c| _ 3fcc 2 , 



Let us consider the change of variables a e = £f sin 2 ij, where £f = a^, r\ — rj (t). Then, because 

,-6/e = Q g 2 M (V + 1 + 

^ " 8^G K + l| 3 / e ^ + 1+ ^J' 



a^ 2 /(3u> + l) / fee 2 \ 2fc, 



a =e lA sin 2 / e 7 yo , 
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Eq. (jlip after simplifications results in 



sin 2 770 = 1/N 



Wp + 1 

w 



2 a/3 



Wq 



3/2£ 



(cos^) 1+3 ^ e (sin 77) 1+3 / e drj 



(Ha V u 'o + 1 J J I ^ ^ gm 6/e ~ + ~ ^ W( ^ ^ e ^ cos 6/ £ rj + i (w Q , a Q , e) sin 2 / e r\ cos 6 / e 77 

where we have introduced the following three parametric coefficients 



dt, (12) 



At(wo,ao,e) 



Wq 



\wq + 1| 3 A 



2fc + 1 



H 2 n 2 
n a 



3w + 1 

Wq 



H a 



2ko 

Wq ' 



l(wo,a ,e) 



ike 2 

HqQq 



Wq 



Wq 



i/'- 



Thus, integration of Eq. (12) from r\ = to r\ = t\q = arcsin ^a/I + ^/wqj gives the age of the universe. The 

lifetime of the universe can be found by integrating Eq. (12) from rj — to rj = Tt/2. Wc shall in what follows consider 
various possible cases. 

1. Flat spacctimc. The coefficients /1 and v in this case depend only on wq and e. Coefficient 7 = 0. The numerical 
calculations of the age of the universe and the difference between time of hnal singularity t f and to for various values 
of wq and e are given in the following table. The used Hubble parameter for that calculation is Hq — 72 km/s/Mpc. 





w = -1.01 


wq = -1.05 


Wq = - 


1.1 


e 


tf -t , Gyr 


h, Gyr 


tf -t , Gyr 


h, Gyr 


tf -t , Gyr 


to, Gyr 


0.1 


252.116 


14.443 


103.026 


13.842 


61.938 


13.245 


0.2 


164.752 


14.442 


76.156 


13.838 


48.557 


13.838 


0.5 


86.958 


14.440 


45.763 


13.827 


45.763 


13.827 


1 


50.798 


14.436 


28.691 


13.813 


28.691 


13.813 


3 


19.726 


14.429 


11.916 


13.771 


11.916 


13.771 


6 


10.286 


14.423 


6.353 


13.756 


6.353 


13.756 


9 


6.929 


14.420 


4.319 


13.745 


4.319 


13.745 



2. Spacetime with positive spatial curvature. The coefficients v, 7 depend on wq, e and ao- We use the values 
1, 5 and 10 for the relation a 2 H 2 /c 2 . We obtain in this case 
a )a 2 H 2 /c 2 = 1. 





w = -1.01 


W q = -1.05 


Wq = - 


1.1 


e 


tf -t , Gyr 


to, Gyr 


tf - to, Gyr 


to, Gyr 


tf -t , Gyr 


to, Gyr 


0.5 


74.752 


11.335 


39.815 


10.941 


27.319 


10.535 


1 


44.344 


11.333 


25.400 


10.931 


18.107 


10.518 


3 


18.017 


11.327 


11.000 


10.905 


8.105 


10.471 


6 


9.734 


11.323 


6.025 


10.885 


4.476 


10.436 


9 


6.671 


11.300 


4.145 


10.874 


3.089 


10.420 



b)a 2 H 2 /c 2 = 5. 
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w = -1.01 


w = -1.05 


Wo — 


i i 

J L 


6 


L f L 0, 1 


b o, ^y L 


ft — tn Gvr 


°0, ^j 1 


tf — to, 'jyr 


+ City- 


0.5 


83.422 


12.335 


43.857 


11.932 


29.912 


11.518 


1 


48.680 


12.333 


27.454 


11.923 


19.471 


11.501 


3 


18.825 


12.327 


11.341 


11.899 


8.359 


11.459 


6 


9.763 


12.324 


6.013 


11.883 


4.491 


11.429 


9 


6.550 


12.288 


4.073 


11.873 


3.060 


11.413 



c)a 2 H 2 /c 2 = 10. 





w = -1.01 


w = -1.05 


w = - 


1.1 


e 


tf - to, Gyr 


h, Gyr 


tf -t , Gyr 


to, Gyr 


tf -t , Gyr 


to, Gyr 


0.5 


85.097 


13.177 


44.749 


12.700 


30.527 


12.216 


1 


49.668 


13.174 


28.022 


12.689 


19.881 


12.197 


3 


19.227 


13.168 


11.594 


12.622 


8.550 


12.148 


6 


9.989 


13.164 


6.159 


12.643 


4.603 


12.114 


9 


6.711 


13.125 


4.177 


12.620 


3.139 


12.086 



3. Spacetime with negative spatial curvature. The coefficients fi, v, 7 depend on wq, e and ao- One notes that the 
right hand side of Eq. (fTTj) is negative for a^Ho/c 2 = 1 along the interval —1.1 < wo < —1.01 because we consider 
the values 5 and 10 forrelation d^H^fc 1 We get 

a)alH 2 /c 2 ^5. 





w = -1.01 


Wo = -1.05 


w = - 


1.1 


e 


tf -t , Gyr 


to, Gyr 


tf -t , Gyr 


to, Gyr 


t f -t , Gyr 


to, Gyr 


0.5 


90.991 


18.850 


47.949 


17.462 


32.756 


16.261 


1 


53.224 


18.842 


30.118 


17.435 


21.405 


16.219 


3 


20.775 


18.826 


12.592 


17.373 


9.305 


16.119 


6 


10.907 


18.817 


6.762 


17.334 


5.060 


16.054 


9 


7.385 


18.800 


4.617 


17.241 


3.471 


15.943 



b)alH 2 /c 2 = 10. 





w = -1.01 


w = -1.05 


w = - 


1.1 


e 


tf -t , Gyr 


to, Gyr 


tf -t , Gyr 


to, Gyr 


- to, Gyr 


to, Gyr 


0.5 


88.946 


16.241 


46.853 


15.381 


31.997 


14.564 


1 


52.009 


16.237 


29.412 


15.362 


20.893 


14.533 


3 


20.267 


16.226 


12.267 


15.317 


9.059 


14.457 


6 


10.613 


16.219 


6.568 


15.287 


4.913 


14.406 


9 


7.170 


16.203 


4.477 


15.212 


3.364 


14.316 



It can be checked that the above numerical estimates provide realistic results for an age of the universe which 
weekly depends on e and wo- The lifetime of universe, which on the contrary strongly depends on e and wo, decreases 
with growing e and Wq. 



IV. NEW SOLUTIONS WITH BIG FREEZE SINGULARITY 

We shall consider in this section a flat universe which shows a big freeze singularity (BFS) in its future, i.e. 

and 

P =-c 2 (4/r 6/e (a}-a £ )" 1_3/e . 

Putting e = 3, we have 

2 a2 3 2 
p = —c p djp . 
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In this case, one can solve the above equations to find the scale factor. We obtain it in parametric form 

I • \ 2 /3 

a = af [smr]) , 



t = t f + ~ (In | tan §| +COS77) , v ; 



where 




cos 77 

at = — : drj. 

k sin 77 

We set now tf = and therefore 77 = corresponds to t = —00, r\ = ir/2 to t = (BFS) and 77 = 7T to t = +00. 
Now one can see that if tp = a 3 then the following equation holds 

^£ = («(*) -AM (14) 

in which 

2k 2 sin 2 77(2 cos 2 77 + 1) , 

w = ■: , A = — 4k . 

cos 4 7/ 

This is very interesting point for the spectral theory of the Schrodinger equation. The potential v(t) — > +0 at 
f — *■ ±00 and v(0) = +00. But wc have a bounded state (tp € L 2 and no zeros at t € (—00; +00)), and this is the case 
notwithstanding for which the potential has a singularity at t = {rj = 7r/2). One can check that 

/oo 
V> 2 cft = 1, 
-00 

if tp — ^/15k/4 sin 2 ry. 

Now we can use Eq. (fT4|) to find the second solution tp with the same potential v and the same value of the spectral 
parameter A, i.e. 



^ = ^ I "I- 



dt 

~2 



We get 



t „ 4 cos 77 „ . 9 . 

= — 2 cos 77 H — 2 sm 77 m 

sin 77 



cot- 
2 



(15) 



Therefore we have a new solution for the same expression 



A= 12ttG[ P - ^ 



This solution describes two universes: the first one begins at t = — 00 and then progressively contracts until a big 
crunch singularity at t = (or 77 = tt/2). The second solution begins at t = (big bang) and then starts expanding. 
One can see that the Hubble root H = dlnip 1 / 3 /dt has the asymptotic behavior given by 

lim H = ±2k, 
f ±00 

and iJ(0) = 00. Thus we have dS universe at t — > ±00. 

A most interesting solution is the superposition of tp = sin 2 77 and ip. We can see that \Jj = c\ip + C2?/> results in a 
new solution a genera i = 'J 1 / 3 , such that 

— = (v(t)-\)V, (16) 

for the same v and A. This solution describes three distinct kinds of universes. If c\ = 1 and C2 = —0.1, then we get 
Universe I. This universe begins at t = — 00 and then progressively contracts until a big crunch singularity at t = tf. 
—00 <ti<0 (for the case c\ = 1, C2 = —0.1 we have ij ~ 0.72). 
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Universe II. This universe begins at t = ti (big bang at a genera i(ti) = 0) after which it starts expanding until a BFS 
which takes place at t = 0. It can be shown that this type of universe cannot be fitted to the data of astronomical 
observations and that it does not lead to a phase of accelerating expansion. 

Universe III. This universe begins at t = as a BFS (with a genera i ^ 0) and then starts expanding until it finally 
behaves like De Sitter for large t 

V M - ?! 

t— >oo o 

We find the latter two kinds of cosmological evolution to be most interesting. 



V. SOLUTIONS WITH COSMOLOGICAL CONSTANT 



Let's consider the universe filled with generalized Chaplygin gas and nonzero vacuum energy. The equation for the 
scale factor 

a 2 8nG. /T 6 / e 

* = — ( (a;- fl ')3/« +A) (17) 

can be conveniently rewritten in terms of the rescaled variables A — a/af and T = (87rG|A|/3) 1 / 2 t. If, without loss 
of generality we set (3~ 6 / e aJ 3 |A| -1 / 2 = 1, we fanally derive from (jXTJ) an equation for A, 

A'" 1 |A| 

A 2 (1 - A e ) 3 / e A ' 

with the prime denoting derivation with respect to T. The solution for this equation can be now given in parametric 
form: 

a) For a universe starting at a big bang 

A = (sin?7) 2 / e , 

T = Tf — 2e~ 1 f* /2 d7/cot77(cos77) 3 / e (l ± (cos r,) 6/e )~ 1/2 , ^ 

where the signs "+" and "-" correspond to positive and negative cosmological constant, respectively. We put Tf = 
and therefore ?y = corresponds to T = — oo (big band) and r\ = tt/2 to T = (BFS). 

b) For a universe starting at a BFS 



A = (sin??) 2 / 6 , 

Setting Tf = and therefore rj = ir/2 corresponds to T = (BFS) and 77 = it to T = —00 (big chrunch) 



T = T f + 26- 1 dr7Cot77(cos77) 3 / e (l ± (cos v) 6/e )' 1/2 , ^ 



VI. SUMMARY AND CONCLUSION 



This paper deals with several new cosmic solutions, which show or do not show a future big freeze singularity, 
derived from a universe filled with generalized Chaplygin in the cases where: (i) dark energy is the sole component of 
the universal energy, (ii) the universe contains in addition some amount of dark matter, and (iii) the universe is also 
equipped with a cosmological constant. We discuss the observational feasibility of such theoretical solutions, reaching 
the conclusion that some nonzero amount of dark matter should be present in order for the given model with future 
big freeze singularity to be compatible with current observations. 

It is worth noticing that the above conclusion and the alluded results are only valid in a classical framework. It 
could well be expected that our classical treatment would break down in the neighborhood of the singularity where the 
energy density tends to infinity. Actually, a quantum-gravity consideration which contemplated the Planck length as 
the ultimate resolution limit would smooth out such a singularity and allowed for a further evolution for the universe. 
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